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A simple lattice gas model in one dimension is constructed in which each site can be occupied 
by at most one particle of any one of D species. Particles interact with a randomly drawn nearest 
neighbor interaction. This model is capable of reproducing the factorial moments observed in high- 
energy scattering. In the limit D oo, the factorial moments of the negative binomial distribution 
are obtained naturally. 

PACS numbers: 13.85.Hd,13.65.+i,24.60.Lz,25.75.+r 
, Factorial moments have provided a useful tool for the analysis of high-energy scattering data as obtained in e+e" 
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scattering , pp scattering (at energies up to 900 GeV) |^ or the scattering of protons or heavy ions by heavy nuclei 
at a projectile energy of 200 GeV/A One considers a range of some variable (usually the rapidity) for which 
one knows the average multiplicity, (n), and its dispersion, (An^). The data is broken into M equal bins, and one 
constructs the factorial moments as 
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^ ■ Empirical factorial moments found in such high-energy scattering experiments are well reproduced by the 'negative 
^sD . binomial distribution' (NB) for which 
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The form of Eq. (|^) invites the consideration of 'universal' plots of Fq versus F2 which permit the comparison of data 
" from very different processes Such plots have been made over the available range 1 < ^2 < 1-8. They reveal 

I both a remarkably universal behaviour and striking agreement with the negative binomial distribution. Giovannini 
and Van Hove Q suggested that these results could be understood 'in terms of partial stimulated emission of bosons, 
or of a simple form of cascade process, or (more artificially) with both mechanisms'. 

The most interesting feature of Eq. (^) and the data which it fits is the growth of the factorial moments with 
^ increasing cM. This has sometimes been regarded as evidence for the presence of fluctuations on many different 
" " ' scales with the related possibility that this may indicate the presence of critical phenomena. Chau and Huang j|] 
have offered an alternate view. They identified the full range of rapidity with the iV sites of a one-dimensional Ising 
(or lattice gas) model. They constructed the factorial moments in this model analytically (for N 00) and set the 
two parameters of the Hamiltonian by fixing the global values of (n) and (An^). The resulting factorial moments 
agree through 0{M) with Eq. (^ and underestimate higher terms. Agreement in the constant and linear term of 
Fq{M) is actually a trivial consequence of (i) the fixing of (n) and (An^) and (ii) the fact that the Ising model predicts 
many-body correlations of a finite range |^ . 

Our purpose here is to demonstrate the surprising result that an elementary extension of the lattice gas model 
permits exact replication of Fq^{M) for all q and M . This extension is simply stated: Construct a one dimensional 
lattice gas model (in the N ^ 00 limit) where each site can either be empty or occupied by one particle which can be of 
any of D species. Each species has a chemical potential, fid, and each pair of species has a nearest neighbor interaction 
of strength e^d'- This model can be solved analytically for the q-body correlation functions and factorial moments 
using textbook techniques. As usual, this involves the construction and diagonalization of a matrix, AA, related to 
the partition function for one pair of adjacent sites. The parameters in this model can be chosen to reproduce the 
results of Chau and Huang for any D. Another choice of parameters will be shown to reproduce the NB results of 
Eq. (|). 

Solution of this model involves a (I? + 1) -dimensional real, symmetric matrix, AA: 
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A^oo = 1 

Mod = cxp [-Md/2] (4) 
Mdd' = cxp [-edd' - A*<i/2 - A*d'/2] 

where matrix indices run from to _D ]To| . We also require the number operator at site i, rii = 1 — T, with 
Tdd' = Sdo^d'o- We shall arbitrarily set all chemical potentials equal to /io- In the thermodynamic limit (both N 
and — !■ cx)), it is sufficient to consider the Z?-dimensional submatrix, M, obtained by neglecting the 0-th row and 
column of M. The coupling of M to the remaining elements of M can then be treated exactly using first order 
perturbation theory. Correlation functions and factorial moments can be determined in terms of the diagonal form, 
Md, and the related orthogonal matrix 6. (Where M is given as OMdO^ ■) For example, the two-body correlation 
function for microscopic sites i and i + j is 
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where the Xk are the eigenvalues of M and the are normalized coefficients which follow from the eigenvectors 
of M as 
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where J\f is chosen such that — 1- The second factorial moment is then given as 

F^^^{M) + (7) 
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After considerable algebraic effort, one finds the following expression for the factorial moments of the generalized 
lattice gas model: 
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Here, we have introduced the definition 
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which implies si — c/2. 

Eqs. (||) and ( p^ ) represent the primary results of our generalized lattice gas model. There are two special cases of 
interest. First, if 
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one immediately obtains the results of Chau and Huang for any D. This can be realized either when the sums of 
Eq. contain only one term or when all the eigenvalues are equal. The second special case corresponds to the 
choice 
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With this choice, Eqs. (|) and (|l^ reduce to the NB resuhs of Eq. (|) 
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It remains to be seen if the parameters in Ai can be selected so that the constraints of Eq. ( p^ are satisfied. The 
following simple prescription works: Set all the off-diagonal elements of A4 equal to zero. Make a draw of D random 
numbers, Xd, from the interval [0,1]. Set Mdd equal to XdL. For each draw, choose L such that si = c/2. (This 
sets the dispersion to its empirical value for each draw.) The physical content of this prescription is clear. Identical 
particles experience a nearest neighbor interaction which ranges from « —fJ^o to +00. Inequivalent particles experience 
a nearest neighbor interaction in the range —fio ^ f-dd' < +00. More 'democratic' schemes can also be constructed. 
They share the feature that M is sparse with roughly D (randomly selected) elements non-zero ||l^ . 

Our prescription meets the remaining conditions of Eq. (|l2|) with increasing accuracy as -D — *■ 00. We shall illustrate 
this with a numerical study of the first seven factorial moments. We take (n) — 20 and (An^) — 110 as would be 
appropriate for the description of pp scattering at 200 GeV. For these data, c{n) = 4.5. (Qualitatively similar results 
are obtained for 900 GeV pp scattering.) We consider the ratios 
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These ratios should be = 1 to reproduce the negative binomial distribution. Our constraint ensures ri — 1. For 
each value of D, we have drawn 10^ matrices according to the prescription above. In Table || we report the 'ensemble 
average over theories', {{rq)), and its dispersion for 2 < g < 6 |jl^ as obtained for D = 1, 2, 4, . . . , 512. Since there 
are no parameters to adjust, this is an extremely stringent test of our prescription. It succeeds. 

Several comments arc in order. For fixed q, the value of {{rq}) approaches 1 like 1/D as D 00. The dispersion 
also vanishes (like I/a/D). Thus, as D becomes large, our simple prescription converges to the results of the NB 
for any fixed q. For fixed D, the error in {{rq)) and its dispersion grow as (7 — > 00. The value of Z) = 16 results in 
sufficiently small errors and dispersions that randomly drawn dynamics have a high probability of reproducing the 
empirical factorial moments for pp scattering at 200 GeV within existing experimental uncertainties. This value of 
Z) = 16 is also sufficient to provide a quantitative description of factorial moments for 900 GeV pp scattering and, 
indeed, of all other high-energy scattering experiments for which factorial moments are known. Our point here is 
that there exists at least one simple prescription for satisfying Eq. (|l2|). Other more efficient prescriptions may well 
exist. 

We have shown that a simple one-dimensional lattice gas model with D species and randomly drawn nearest neigh- 
bor interactions between equivalent species can reproduce the factorial moments of the negative binomial distribution 
as D becomes large. Given the limited range of cM covered by current experimental data, a remarkably small number 
of species is sufficient to provide a quantitative description of the empirical factorial moments. This offers some 
understanding for both the success of cascade calculations and the anecdotal observation that the results of such 
calculations are often surprisingly insensitive to the details of the model. In the present picture, any of our randomly 
drawn theories would also be likely to succeed (at least at the level of the factorial moments) . 

The empirical observation that factorial moments grow like powers of cM in pp, e^e~ and relativistic heavy ion 
collisions (the phenomenon of intermittency) has often be taken as evidence of the existence of fluctuations on 'all 
length scales'. As such, it is sometimes seen as an indicator of the presence of a non-equilibrium, critical phenomenon. 
While we do not deny the possibility that intermittency can be a signature of critical phenomena, we have shown 
that a simple but highly heterogeneous (equilibrium) system can also lead to intermittency. In short, intermittency 
is not a unique signature of critical phenomena. Given the small number of species required by our model to fit the 
factorial moments obtained in high-energy pp scattering, a critical phenomenon docs not even appear to be the 'most 
plausible' cause of intermittency. 

Two of us (AD J and TW) would like to acknowledge the hospitality of NORDITA. This work was partially supported 
by the U.S. Department of Energy under grant no. DE-FG02-88ER 40388. 
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TABLE L The ensemble average and dispersion of as defined in Eq. ( [l3| ) for 2 < g < 6 and various values of the number 
of particle species, D. The average was taken over 10^ randomly drawn theories with c{n) = 4.5 for each draw. The case D = 1 
corresponds to the ordinary lattice gas (Ising) model H and has no dispersion. 
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1 


0.75 


0.5 


0.3125 


0.1875 


0.109375 


2 


0.812 ±0.023 


0.598 ± 0.035 


0.417 ±0.037 


0.280 ± 0.033 


0.184 ± 0.027 


4 


0.881 ± 0.046 


0.723 ± 0.078 


0.569 ±0.091 


0.435 ± 0.091 


0.326 ± 0.083 


8 


0.941 ± 0.067 


0.851 ±0.132 


0.751 ±0.179 


0.652 ±0.207 


0.558 ±0.218 


16 


0.981 ±0.076 


0.951 ±0.173 


0.916 ±0.270 


0.879 ± 0.360 


0.842 ± 0.441 


32 


0.998 ± 0.069 


0.999 ±0.170 


1.005 ±0.295 


1.017 ±0.442 


1.039 ± 0.612 


64 


1.002 ±0.050 


1.007 ±0.128 


1.019 ±0.229 


1.040 ± 0.355 


1.071 ±0.515 


128 


1.001 ±0.035 


1.005 ± 0.088 


1.013 ±0.154 


1.026 ± 0.233 


1.044 ± 0.327 


256 


1.001 ± 0.024 


1.003 ± 0.061 


1.007 ±0.105 


1.013 ±0.155 


1.023 ± 0.210 


512 


1.000 ±0.017 


1.002 ± 0.042 


1.004 ± 0.072 


1.007 ±0.106 


1.012 ± 0.142 
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